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Superconducting vortex lattices, glasses and liquids attract great interest as model systems of 
crystallization and as a source of microscopic information of the nature of superconductivity. We 
report for the first time direct microscopic measurements of the vortex lattice tilt modulus C44 in 
ultra-pure niobium using time-resolved small angle neutron scattering. Besides a general trend to 
faster vortex lattice dynamics for increasing temperatures we observe a dramatic changeover of the 
relaxation process associated with the non-trivial vortex lattice morphology in the intermediate 
mixed state. This changeover is attributed to a Landau-branching of the Shubnikov domains at the 
surface of the sample. Our study represents a showcase for how to access directly vortex lattice 
melting and the formation of vortex matter states for other systems. 

PACS numbers: 71.18.+y, 74.25. Dw, 74.25. Q-(-, 78.70.Nx 



I. INTRODUCTION AND MOTIVATION 

The morphology of superconducting vortex lattices 
(VL) attracts great interest as a source of microscopic 
information of the nature of the superconductivity and 
as model systems of condensed matter. The elastic ma- 
trix f^ap of a VL thereby describes the energy associ- 
ated with a distortion of the VL due to thermal fluctu- 
ations, gradients of magnetic field or temperature, pin- 
ning and the presence of transport currents. Analogous 
to crystal lattices the elastic matrix $c<^ of a VL deter- 
mines the thermal stability and the state of aggregation 
of superconducting vortex matter: Besides the regular 
Abrikosov VL, VL Bragg glasses, liquids and ices have 
been identified^ ^. 

At VL melting transitions the shear modulus discon- 
tinuously jumps to cgg = 0, where the long range order 
vanishes. VL melting was observed for various super- 
conducting systems, mostly for high-Tc compounds due 
to their high transition temperaturd^, but also in com- 
pounds characterized by disorder as NbSe2 and MgBjpM^l 
Surprisingly, a melting transition was recently reported 
also for the heavy fermion compound URu2Si5P in the 
clean limit. Melting transitions show up as characteris- 
tic dips of the differential resistivity and can also lead to 
tiny jumps of the local magnetization of the order of few 
tenths of a Gauss, detectable with sensitive Hall probes^. 
Measurements to detect melting transitions are intricate, 
as effects induced by depinning can yield similar results. 
However, the presence of pinning is required in general, 
as a perfect, pinning free VL shows no signature in the 
resistivity at the melting transitiorP. 



Moreover the elastic matrix of VLs is intimately 
related to the pinning and depinning properties of su- 
perconducting vortices. This is especially important for 
technical applications: If transport currents are applied 
to superconducting materials, the Lorentz force acting on 
the vortices leads — with increasing current — to dissi- 
pative processes such as vortex creep^ , thermally assisted 
fiux-fiow (TAFF^pSl and flux-flow (FFjEl. Therefore, the 
ability of superconducting materials to carry large trans- 
port currents for technical applications is closely con- 
nected to the pinning properties of superconductors and 
the elasticity of the VL. The elastic constants of the VL 
Cii for compression, C44 for tilt and Cge for shear hence 
reflect the microscopic nature of the superconductivity 
as well as impurity or surface proper ties of the supercon- 
ducting sample due to pinnin g^^ttis i 

The experimental access to the elastic matrix <I>q^ 
of VLs — in particular for non-equilibrium states — 
by macroscopic bulk techniques such as the transport 
propertieJi^, the magnetizatiorP or measurements us- 
ing vibrating reedJ^^ is strongly influenced by parasitic 
pinning effects as well as geometric effects. Moreover, 
the momentum-dependence of the elasticity of VLs can- 
not be determined unambiguously by macroscopic mea- 
surements. Due to the use of thin film samples, micro- 
scopic surface sensitive techniques such like decoration 
or magneto-optical methods suffer from similar or even 
stronger pinning and geometry induced effects. In con- 
trast, local probes such as muon spin relaxation (/iSR) 
and scattering techniques such as neutron scattering yield 
microscopic information on bulk VLs. However, as the 
accessible timescale of inelastic neutron scattering tech- 
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niques is still too short for VL dynamical properties, neu- 
tron scattering was up to now limited to characterize the 
static properties of VLs with only a few exception^i^ttSO] 
In this paper we report direct microscopic measure- 
ments of the VL tilt modulus C44 with drastically re- 
duced limitations due to surface pinning in ultra-pure 
bulk Niobium (Nb) using a time-resolved neutron scatter- 
ing technique as combined with a tailored magnetic field 
setup. With its low Ginzburg-Landau parameter k, situ- 
ated close to the border of type-I and type-II behaviour, 
the superconductivity in Nb is ideally suited as model 
system for systematic studies of vortex mattei^^iHSll, By 
imposing a periodic tilting of the magnetic field, we in- 
duce a relaxation process of the VL which can be de- 
scribed by a diffusion process in the limit of uniform tilt. 
The diffusion constant of this diffusion process is given 
by the tilt modulus C44 of the VL and the flux flow resis- 
tivity pff- The characteristic properties of the diffusion 
process are observed by means of time resolved strobo- 
scopic small angle neutron scattering (SANSJP^. The re- 
laxation processes observed show increasing VL stiffness 
with increasing magnetic field H and reduced damping 
with increasing temperature T. This behaviour agrees 
well with calculations performed within a VL diffusion 
modeP^. Besides these general trends, we observe a dra- 
matic change of the relaxation processes associated with 
the non-trivial VL morphology in the intermediate mixed 
state (IMS). 

Our study represents a showcase for how to access di- 
rectly VL melting, the formation of vortex glass states 
and slow vortex dynamics also in unconventional su- 
perconductors, notably the cuprates, heavy-fermion sys- 
tems, borocarbide or ironarsenide systems. 

The outline of this paper is as follows: In the follow- 
ing section (|ll]), we briefly revisit the salient features of 
the elasticity of superconducting VLs. The experimental 
setup which was developed for our measurements is dis- 
cussed in depth in section (III). The results obtained for 
the tilt modulus C44 are presented in section ( |IV[ ). The 
response of the VL to a changed magnetic held environ- 
ment in the k = limit with a diffusion Ansatz^ 10 26 27] 
is discussed in section (|v|. Finally, the relevance of our 
experimental setup for the investigation of different mag- 
netic systems is discussed in section (IVII). 



II. ELASTICITY OF VORTEX LATTICES 

In the following, we introduce the elastic proper- 
ties of VLs. We first concentrate on the derivation 
of the free energy of an arbitrarily curved VL using a 
Ginzburg-Landau ansatz. We then derive the elastic ma- 
trix and deduce the eigenfrequencies, characteris- 
tic timescales and the moduli for compression, tilt and 
shear of superconducting VLs. In particular, we discuss 
the elastic moduli for uniform tilt which are relevant for 
the interpretation of our measurements of the VL elas- 
ticity in superconducting Nb. For a detailed theoretical 
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FIG. 1: Panel (i) schematically depicts the inter (dotted) and 
intra (dashed) vortex interactions determining a particular 
symmetry, structure and elasticity of VLs. Panel (ii) depicts 
the cusp-like response of a superconducting VL to a pinning 
force dF. Panel (iii) schematically depicts a thermally fluctu- 
ating superconducting VL for a temperature below the melt- 
ing temperature Tm- 



description on VL elasticity, we refer tcPiSEMZl, 



A. Elastic Energy of Vortex Lattices 

Similar to the elasticity of crystal lattices, which is de- 
termined by electrostatic or covalent forces, the elastic- 
ity of superconducting VLs is determined by the vortex- 
vortex interactions. The free energy F of an arrangement 
of arbitrarily curved vortices in the Ginzburg-Landau 
regime close to Tc has been approximated by BrandiP 
and can be written in terms of two contributions 
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with 



A' = A/(|^P)V2^A/(1-6)V2, 

e' = e/[2(i-&)]^/^ 



(1) 



(2) 



(3) 



(4) 



and the reduced field b — B/Bc2- The first term of the 
sum in eq. ([T]) represents the repulsive electromagnetic 
vortex-vortex interaction with an effective London pene- 
tration depth A'. The second term represents the attrac- 
tive interaction of vortex lines due to the condensation 
energy of overlapping vortices with an effective coher- 
ence length The vortex self energy or line tension is 
included in the diagonal terms i = j. A schematic sketch 
of the interactions is given in Fig. [l] panel (i). 
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FIG. 2: (i) Schematic depiction of the elastic constants cn for 
compression, C44 for tilt and cee for shear of a superconducting 
VL. (ii) Depiction of the k dependence of the elastic matrix 
^a,p and of the elastic constants cn, C44 and cge of the VL. 
Note the different scaling for ceJ^. 



The elastic energy of a distorted VL, caused by pin- 
ning, structural defects, field gradients or transport cur- 
rents, temperature gradients or thermal fluctuations, is 
small for most cases. Therefore, it can be calculated 
by linear elastic theory expressed in fc-space. The dis- 
placements of a vortex line Ui{z) = rj(z) — 'R.i{z) — 
(ui^x'iUi^y^O) from its ideal position = {Xi;Yi;z) is 
expressed by its Fourier components 
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With u(k) = (u^; Uy] 0) the elastic free energy reads 



2 Jbz 8^ 



ia(k)$a,3(k)M«(k) 
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where {a, (5) = {x,y). The integrals in eq. ^ and eq. 
^ cover the first Brillouin zone of the VL in k-space 
and <kz< C^^, respectively. $Q^(k) is called the 

elastic matrix of the VL. 

$Q./3(k) is real, symmetric and periodic in A:-space and 
is related to the elastic moduli cn for compression, C44 
for tilt and Cgg for shear within continuum theory by 



$Q/}(k) = (cu - C(iQ)kakf}+ 

with k\ = {kl + fcp (cf 
of the elastic matrix is plotte 
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dependence 
panel (ii). 



The Labusch parameter describes the elastic inter- 
action of the VL with pinning potentials caused by ma- 
terial inho mog eneities . For individual pinning, is k- 
independentPSl, for weak collective pinnin^^Si^ aL(k) de- 
creases when fcj^ > or k^ > where Rc and 
Lc = {c4i/ceey^^Rc are the radius and length of the 
coherent short range ordered regions of the pinned VL. 



B. Uniform Distortions 

For uniform distortions, the elastic moduli of the VL 
may be written ad3 
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with the Ginzburg-Landau parameter k. Cn — cge is the 
modulus for isotropic compression. H is the applied field, 
which is in equilibrium with the VL at the equilibrium 
induction B, given by the magnetization M = ^qH—B < 
0. The response of a VL to a change of the magnetic 
field direction is characterized by the tilt modulus C44 
of the VL. Note that cge vanishes either for B — Bc2, 
which corresponds to strongly overlapping vortex cores, 
for A — ^ 00, corresponding to strongly overlapping vortex 
fields, or for k = l/v^- In the special case k — l/\/2, all 
VL arrangements have the same free energj^^. 

Due to the long effective interaction lengths A' and 
Cn and C44 strongly depend on the /c- vector of the 
disturbance which is referred to as non-locality of the 
VL. This leads to a strong softening of the VL for short 
range distortions. This non-locality gives rise to large 
distortions, caused by pinning, disorder or thermal fluc- 
tuations, whereby the VL reacts to external forces in the 
form of a sharp cusp and not like pulling a string (cf. 
Fig.[l]). 



C. Characteristic Timescales 

The elastic matrix of superconducting VLs determines 
the restoring force of k-dependent distortions. The eigen- 
modes and the characteristic damping of VL fluctuations 
is determined by the restoring force as well as the vis- 
cous damping of the VL motion. The movement of the 
vortices with the velocity v is damped by the viscosity 
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This creates a drag force vi] per unit volume, where ppF 
represents the flux-flow and p„ the normal conducting re- 
sistivity. The elastic eigenmodes of the VL are given by 
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a diagonalization of the elastic matrix $0/3. The result 
is a compressional and a shear eigenmode, relaxing with 
exponential time dependencies. In continuum approxi- 
mation, this yields: 

ri(k) = (cn(k)fci + C44(k)fcf )/r; « 
r2(k) = (cgefci + C44(k)fc2)/^ « rifc2/A:2 

For the VL in a typical clean low k superconductor as 
e.g. Nb, the the characteristic relaxation rates Fi and F2 
of the VL are in the range of 10^^ s^^. 

III. EXPERIMENTAL SETUP 

To measure the VL tilt modulus C44 by means of time 
resolved SANS, a time varying magnetic field setup, con- 
sisting of two orthogonal pairs of Helmholtz-coils has 
been designed. By imposing a periodic tilting of the mag- 
netic field, we induce a relaxation process of the VL which 
can be described by a diffusion process. The relaxation 
of the VL is measured by means of time resolved stro- 
boscopic SANS. The measurements have been perfor med 
on the small angle diffractometer V4 at BENSCpl2ni. We 
introduce the details of the experimental setup used for 
our study in the following paragraphs: 

A. Sample Environment and Magnetic Field Setup 

A schematic drawing of the experimental setup is given 
in Fig. [3] The sample is located in the center of two 
Helmholtz-coils, cooled with a closed-cycle cryostat to a 
minimum temperature of 4K. Both magnetic fields and 
the sample can be rocked together with respect to the 
vertical z-axis. The static main field Hgtat, applied along 
the y-axis is generated by bespoke water cooled copper 
coils'^"'^. Hstat is oriented approximately parallel to the in- 
coming neutron beam. The angle enclosed between Hgtat 
and the incoming neutron beam, which is also applied in 
the xy plane, is denoted rocking angle 0. A magnetic 
field 75 mT< Hgtat <135mT was applied. 

The time varying field Hose is generated by a small 
air-cooled set of Helmholtz-coils inside the main coil, 
driven with an arbitrary waveform generator and an am- 
plifier. Hose is oriented along the x-axis perpendicular to 
Hstat- A rectangular pulse shape with an amplitude of 
Hose — mT O 5 mT and a repetition rate of 0.2 Hz was 
applied. The resulting field Htotai is rotated with respect 
to Hstat by the angle e = arctan |Hose|/|Hstat| in the xy 

plane. Hstat > Hose yiclds that |Htotal| « |Hstat| = |H|. 

Due to the perpendicular alignment we omit the vectorial 
notation of Hstat and Hose- A smearing of the applied 
pulses is caused by the rise and fall time of amplifier used 
for iJose- It has been determined with a Hall- probe at the 
sample position and found to be ^ 5 ms. 

Hence two equilibrium positions for the magnetic field 
and also the VL emerge, which are separated by e « 2° if 
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FIG. 3: Schematic depiction of the experimental setup used 
for time resolved SANS measurements of the VL tilt modulus 
C44. The sample is located in the centre of two orthogonal 
magnetic fields, generated by Helmholtz coils. With a combi- 
nation of a static magnetic field //atat || y and a time-varying 
magnetic field Hose || x with Hose ^ Hstat, the resulting mag- 
netic field can be rotated with respect to the sample in the 
xy plane. Both magnetic fields and the sample can be rocked 
around the vertical 2;-axis with respect to the incoming neu- 
tron beam by the angle </>. The resulting magnetic field Htotai 
is roughly parallel to the incoming neutron beam. The scat- 
tered intensity is recorded on a two-dimensional detector. 



Hose is alternated between Hose —0 mT and i?ose =5 mT. 
The cylindric sample is aligned with its symmetry axis 
parallel to the z-axis, i.e. a constant demagnetizing factor 
N = 1/2 applies for all angles e. 

For the magnetic field range of our experiment the VL 
assumes a six-fold scattering pattern with a Bragg angle 
of a few tenths of a degree. The instrumental resolution 
in the direction perpendicular to the scattering vector q 
yields a value of 0.2° for a coUimation length of Li = 8 m, 
a sample detector distance of L2 = 8 m, a source aperture 
of i?i = 10 mm and a sample aperture of R2 = 2 mm, 
respectiveljSU. The shift of the VL direction of e « 2° 
is hence much larger than the instrumental resolution 
giving rise to a large contrast. 

B. Measurement Principle 

In the following we describe the principle used for the 
measurement of the VL motion and relaxation, driven 
by the time dependent transverse magnetic field. The 
rocking angle (f) is initially adjusted to satisfy the Bragg 
condition for a reciprocal lattice vector of the VL q = 
GvL = ki — k/ (lying in the xy-plane) for Hose — OmT, 
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i.e. e = 0°. The observed scattering intensity at the 
2D detector at Gvl thus is a measure for the quantity 
of VL which points in this direction. i?osc oscillates be- 
tween /iQ-ffosc = niT and /io^^osc = 5 niT. Thus, the re- 
laxation process between these two equilibrium positions 
can be followed by measuring the integrated intensity at 
the Bragg reflection at Gvl as function of time. Two 
different time-dependent processes can be measured. 

• The first process (^o-f^osc decreases from 5mT to 
OmT) is attributed to the VL relaxing into the 
Bragg conditiorP^. The corresponding time-scale 
is denoted ri. 

• The latter process (/xo^^osc increases from mT to 
5 mT) describes a time-scale, necessary to pull the 
VL out of the Bragg condition denoted as T3. 

In the following, this method is denoted fixed angle 
scan. Typical data for a temperature of T = 4K and a 
magnetic field of fioH =100 mT is shown in Fig. |4j panel 
(i), where the red line indicates the modulus of Hose- 
The whole relaxation process of the VL can be traced 
angle and time resolved, when such scans are performed 
for each rocking angle 0. The latter method is denoted 
time resolved mapping. A representative scan is given in 
Fig.|4) panel (ii) for T = 4K and ^iqH = lOOmT. Fixed 
angle scans are represented by cuts at a fixed rocking 
angle (f) in the time resolved mappings. It is important 
to note that the angular distribution of the VL is always 
integrated over the complete sample and thus is addition- 
ally convoluted with the intrinsic VL mosaicity and the 
angular resolution of the small angle scattering instru- 
ment. 
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FIG. 4: Panel (i) shows a typical fixed angle scan. A trans- 
verse magnetic field defiects the direction of the VL in the 
sample leading to a drop in the Bragg intensity. The red line 
indicates the modulus of -ffosc ■ Panel (ii) shows a typical time 
resolved map for varying rocking angles tj). Two equilibrium 
positions for the VL, as induced by the change of Hose, are 
visible for (j> — 0.5° and (j) = 2.5°. The switching process 
between these equilibrium positions can be monitored as a 
function of time. The horizontal white line in panel (ii) rep- 
resents the fixed angle scan given in panel (i). Both, scans (i) 
and (ii) have been performed at T = 4K and /xo-ff =100 mT. 



C. Stroboscopic Small Angle Neutron Scattering 



The fundamental principle of stroboscopic neutron 
scattering is the excitation of the sample by an exter- 
nal control parameter followed by a measurement of the 
time dependent dynamic response and relaxatiorP^, in 
our case the direction of the magnetic field. To increase 
time resolution and signal statistics, these measurements 
are performed in a stroboscopic manner, i.e. the measure- 
ment is repeated many times where the data obtained for 
the individual cycles is summed coherentlji^H. The stro- 
boscopic SANS technique is realized, using a standard 
SANS setup, extended by a time resolved position sensi- 
tive detector. The repetition cycles of the time resolved 
detector and the control parameter are phase locked. 

The time resolution is mostly determined by a smear- 
ing of the single frames, caused by the wavelength spread 
A A/ A of the neutron beam leading to a variation of 
the neutron time of flight from the sample to the de- 
tector. The time of flight is given by the equatiorP^ 
tTOF[nis] = A[A] • L2[m] • 0.253. A represents the wave- 
length of the neutron beam and L2 the sample-detector 
distance, respectively. For SANS measurements, L2 and 



A determine the accessible g-range. A large sample- 
detector distance L2 and large wavelength A is desir- 
able to resolve large real space structures associated with 
small g- vectors. This leads to a significant loss in time 
resolution. For the measurements on the VL in Nb, 
presented in this manuscript, a wavelength A =8A , a 
wavelength spread AA/A = 0.1 and a detector distance 
L2=8m lead to a time resolution of ±1.6 ms. 



D. Experimental Parameters 

Measurements of the tilt modulus C44 of the VL have 
been performed for an applied magnetic field I-IqH = 
75 mT, 100 mT and 135 mT, each for sample tempera- 
tures between 4K and T^. A schematic sketch of the 
measurement range with respect to the phase diagram, 
obtained in previous SANS measurement^^!! jg given in 
Fig. [5j panel (i): Both, the IMS and the crossover to the 
Shubnikov phase are covered, where the vortex-vortex 
interaction changes from attractive to repulsive. 
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FIG. 5: Panel (i) depicts the phase diagram obtained by previ- 
ous SANS measurement^^. A demagnetizing factor N = 1/2 
applies. The dashed red lines locate scans, performed for our 
measurements of C44, in the phase diagram. Panel (ii) depicts 
the intrinsic mosaicity of the VL as measured with SANS as 
function of temperature for magnetic fields ^qH = 75 mT, 
100 mT and 135 mT. The lines serve as guide to the eye. 



The intrinsic mosaicity of the VL as inferred from 
static rocking scans is given in panel (ii) of Figjsj For 
the highest temperatures the mosaicity of the VL is lim- 
ited by the instrumental resolution. For decreasing tem- 
peratures T < 6K for ^iqH = 75mT and T < 5K 
for ^.qH = 100 mT the mosaic spread shows a pro- 
nounced increase. However, in the Shubnikov phase for 
^{)H =135 mT only a weak increase of the mosaic spread 
is observed for decreasing temperatures. The increase of 
mosaicity for decreasing temperature is attributed to the 
crossover to the IMS where an additional bending of vor- 
tices is caused by the complicated VL domain structure 
in combination with increasing Meissner effect. 

To avoid hysteretic effects, all measurements have been 
taken after cooling in a field to the desired temperature 
(FC). Furthermore, the direction of the magnetic field 
is oscillating continuously due to i?osc- This leads to 
an effective depinning of the VL. An equilibrium state 
can thus be assumed. In analogy to HFC paths in the 
B-T phase diagram of type-I superconductors in the in- 
termediate state, where the magnetic flux is expelled for 
decreasing field, a similar behaviour and morphology is 
expected in the IMS for type-II superconductors upon 
FC. The result is supposed to be an open, multiply con- 
nected topology of Shubnikov domains enclosing regions 
of Meissner phase. In addition, the IMS is characterized 
by Landau branching of the Shubnikov-domains at the 
surface of the sample. 



E. Sample Used for our Study 

For our studies, a cylindrical Nb single crystal with a 
length of 20 mm and a diameter of 4.5 mm was cut by 
spark erosion from a rod that had been produced at the 
ZFW Dresden more than 30 years ago. The cylindric 
symmetry axis of the sample coincides with a crystallo- 
graphic (110) axis. A further (110) axis is oriented par- 
allel to the incident neutron beam. The preparation pro- 
cess of the rod consisted of purification by liquid-liquid 
extraction combined with chlorination and thermal de- 
composition of NbCls followed by electron beam floating 
zone melting, decarburization in oxygen atmosphere and 
annealing in UHV'^'^'^^. The impurity content was esti- 
mated to be less than 1 ppm for interstitial and better 
than 2 ppm for substitutional impurities. 

To remove the surface layer of the cutted sample and 
to decrease the surface roughness, the sample was etched 
with a mixture of concentrated HF and HNO3 for sev- 
eral minutes. In order to remove interstitials, in particu- 
lar hydrogen introduced due to the spark erosion cutting 
process as well as the long storage, the sample was again 
RF-annealed in UHV above 2000° at the University of 
Birmingham for one week followed by surface oxygena- 
tion to reduce the Bean-Livingston barrier for surface 
pinninJ^H. The sample shows a smooth, highly polished 
surface with a pale golden colour without any signs of 
pores, scratches or damages. 

The residual resistivity ratio (RRR) of the sample was 
measured with an eddy current decay method at the Uni- 
versity of Birmingham. The RRR was extrapolated to 
T — > at i? = for a temperature range of 14.5 K 
- 9.3 K assuming Pphonon oc as well as extrapolated 
to B = 0(T = 4.2 K) from Bc2(T = 4.2 K) assuming 
Pb oc yielding values from RRR=8000 to 16000, re- 
spectively. However, the temperature extrapolation is 
more reliable, leading to a value of RRR~ lO''. This leads 
to a Ginzburg-Landau coefficient k ~ 0.74 at 0.9T/T^. 
The AC susceptibility and the magnetization, measured 
at the Technische Universitat Miinchen were consistent 
with the literature. Previous small angle neutron scat- 
tering measurements on the sample have shown no in- 
dications of trapped fiux for decreasing magnetic field, 
further indicating an excellent sample qualitjl^ll and a 
vanishing surface pinning barrier for the vortices. 

The fiux-fiow resistivity pff, responsible for the damp- 
ing of vortex motion in superconductors is related to 
the normal conducting resistivity p„ by eq. (10 1. For 
high purity Nb in the clean limit /3„ was numerically 
approximatecPSl by 



Pn(T) =po + Cn y'T'^ + 



7.212 Je^/j, (e--l)(l-e-) 



dx- 



124.4 ye_/^ (e--l)(l-e-) 



dx 



(12) 
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with b = 1.63 • lO-'^K-^, c = (2.864 ± 0.003) • IQ-^K-^, 
d = (1.81 ± 0.004) • 10-i°K-5 a Debye temperature of 
Qjj — 270 K, Gmin = 35 K and the normahzation con- 
stant c„ = 1.8 • 10-''. The RRR of 10'' and the Ut- 
erature value for p„(300K)~ 1.3 • lO'^f^nPI lead to 
Pa ^ 1.3 • 10^^^ rim. This yields an increase of the nor- 
mal conduction resistivity p„ in the relevant temperature 
range of the experiment from 4 K to 8 K by a factor of 
about 2 from 1.8nJ7cm to 3.2nricm. 



IV. EXPERIMENTAL RESULTS 

A. Time Resolved Mappings 

In the following we present our data obtained for the 
VL relaxation and diffusion. First, we focus on the time 
resolved mappings. Fig. |6] depicts the relaxation of the 
VL at a magnetic field /j.o^^=10GmT for a temperature 
of T = 6.5K (panels (i) and (ii)) and T = 4K (panels 
(iii) and (iv)). The integrated Bragg intensity is plot- 
ted on a linear scale. The horizontal broken white lines 
marked with the black arrows indicate the time when the 
magnetic field direction is switched between the two equi- 
librium positions. The time-range displayed corresponds 
to the gray shading in Fig.|4] panel (ii), however, the 
axes have been rotated for better visibility. Both equi- 
librium positions (indicated with vertical broken white 
lines) at cj) — 0.5° and (p = 2.75° are clearly visible for 
the data obtained at T = 4K. The measurement range 
was reduced for T = 6.5 K due to the limited beamtime. 
Salient features of the VL relaxation are: 

• A larger intrinsic mosaicity of the VL is observed 
for T = 4K in comparison to T = 6.5 K. The addi- 
tional mosaic spread, caused by a possible bending 
of the vortices due to the domain structure of the 
IMS in combination with demagnetizing effects is 
in agreement with static data as displayed in panel 
(ii) of Fig. [5} 

• Relaxation time ti : As expected, the general trend 
is a relaxation characteristic of an exponential de- 
cay (X e-*/^i for both T = 4K and T = 6.5 K, 
indicated by the continuous red lines. The corre- 
sponding time constant ri, obtained for T = 6.5 K 
yields ti = 0.04± 0.008 s, whereas n = 0.16± 0.02 s 
for T = 4K. In addition, we further define r{ as the 
characteristic time-scale, when the VL has reached 
its new equilibrium position. For T — 6.5 K yields 
t[ = 0.2 ± 0.05 s, whereas t[ = 0.45 ± 0.1s for 
T = 4K. 

• Relaxation time T3: T3 indicates the time, at which 
the intensity at the Bragg spot Gvl has decreased 
to 1/e of its initial value. We obtain — 0.06 ± 
0.008s for T = 4K and rg = 0.025 ± 0.005s for 
T = 6.5 K, respectively. 




FIG. 6: Panels (i) and (ii) depict time resolved mappings 
for an applied magnetic field /io-ff=100 mT and a tempera- 
ture of 6.5 K. Panels (iii) and (iv) depict similar scans for 
fj,oH=100mT and 4K. Note the reduced measurement range 
for T — 6.5 K. The contours are plotted on a linear scale. 
The change of magnetic field direction is indicated by hori- 
zontal broken white lines marked with black arrows, whereas 
the equilibrium positions of the VL are marked with vertical 
broken white lines. The continuous red lines indicate the re- 
laxation process of the VL. The red lines serve as guide to the 
eye. 



• The intensity map, obtained for T = 4K is charac- 
terized by a drastic increase of mosaic spread imme- 
diately after the magnetic field direction is changed. 
The time-scale observed for this feature is found to 
be less than T3. In contrast, only a moderate in- 
crease of mosaic spread is observed for T = 6.5 K. 



B. Fixed Angle Scans 

New light is shed on the details of the VL relaxation 
process, focusing on the fixed angle scans as function 
of applied magnetic field H and temperature T. The 
data is given in Fig. [t] panels (i) and (ii) for a mag- 
netic field of fioH = 75 mT, Fig. [sj panels (i) and (ii) 
for pqH = 100 mT and Fig.joj panels (i) and (u) for 
fj-oH — 135 mT. The time-range displayed for the three 
figures corresponds to the gray shadings in Fig. [4j panel 
(i)- 

We first focus on panels (i) of Figs. [7j |8] and [9] which 
show fixed angle scans for increasing sample tempera- 
ture from T = 4K to Tc (as labeled in panels (ii)). The 
increase of scattering intensity after the vertical line, la- 
beled with trigger is attributed to the relaxation of the 
VL into the Bragg condition, the characteristic relax- 
ation time was denoted ri. However, it turns out that 
the qualitative shape of the relaxation process changes 
as a function of temperature. It is therefore not pos- 
sible to determine the characteristic time-scale ri, in a 
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(i) Trigger: H^^_.=5mT to OmT 
Relaxation of VL into Bragg condition 



Trigger: H^^_.=OmT to 5mT (ii) 
Pulling VL out of Bragg condition 



75mT 



7.25K 




'"km 



200 400 600 800 
Time (ms) 



100 200 300 
Time (ms) 



FIG. 7: Panel (i) shows fixed angle scans for hqH =75 mT 
for increasing temperature as labeled in panel (ii). The in- 
crease of scattering intensity after the trigger is attributed to 
the relaxation of the VL into the Bragg condition as func- 
tion of time. The characteristic relaxation time is denoted 
Ti, indicated with black markers. Panel (ii) shows fixed an- 
gle scans for /xo-ff =75 mT for increasing temperature. The 
decrease of intensity after the trigger is attributed to the VL 
being pulled out of the Bragg condition. The characteristic 
relaxation times are denoted T2 and T3, indicated with green 
and blue markers, respectively. Note the different timescale 
for panel (i) compared to panel (ii). The different plots have 
been shifted vertically for clarity reasons, the lines serve as 
guide to the eye. The phase region of the IMS in indicated 
on the rhs. 



way similar to the time resolved mappings^'^ . In order to 
cover the systematic trends, the magnetic field and tem- 
perature dependence of r{ was analyzed. Salient features 
are: 

• The same general trend as observed for the time re- 
solved mappings can also be identified: Increasing 
temperature and increasing magnetic field lead to a 
significantly faster relaxation. The relaxation time 
t[ is indicated with a black marker in each fixed an- 
gle scan. The resulting magnetic field and tempera- 
ture dependence of t{ is shown in Fig. 10 panel (i): 
For low temperatures T « 4 K and magnetic fields 



of noH = 75 mT and 100 mT, t[ =0.5±0.1s. For 
increasing temperatures, t[ decreases characteristic 
of a smooth crossover to values of t[ =0.2±0.05s. 
The crossover temperature thereby decreases with 
increasing field from T = 6.5 K at /xoi/=75mT 
to T = 5.5 K at ^iqH = lOOmT. In contrast, 
for a magnetic field of ^x^H — 135 mT, t[ as- 
sumes a constant value of ri = 0.2 ± 0.05s except 
at the lowest temperature where a slight shift to 
t[ = 0.25 ± 0.05 s is observed. The crossover tem- 
perature for ^iqH — 135 mT is expected to be at 
« 3.5 K, i.e. slightly below the accessible temper- 



ature regime. The magnetic field dependence of 
the crossover temperature of r{ corresponds to the 
crossover from the IMS to the Shubnikov phase. 
The data points obtained for ri from the time re- 
solved mappings are given in Fig. 10 panel (i) as 
well. 

(i) Trigger: H^^^=5mT to OmT Trigger: H^^^=OmT to 5mT (ii) 

Relaxation of VL into Bragg condition Pulling VL out of Bragg condition 
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FIG. 8: Panel (i) shows fixed angle scans for figH =100 mT 
analogous to Fig. [7] 



(i) 



Trigger: H =5mT to OmT 



Trigger: H =OmT to 5mT (ii) 
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FIG. 9: Panel (i) shows fixed angle scans for /lo-ff =135 mT 
analogous to Fig. [7] 



• The characteristic shape of the relaxation changes 
as a function of temperature and magnetic field: 
For the low temperature data points at each mag- 
netic field, namely for 4K < T < 6K at figH — 
75mT, for 4K < T < 5K at fioH = lOOmT and 
for T = 4K at /igi? = 135 mT, the relaxation is 
characterized by a sharp kink, associated with a 
distinct increase of intensity immediately after the 
magnetic field is changed. This sharp increase is 
more pronounced at lowest temperatures. It is as- 
sociated with the step-like broadening of the VL 
mosaic, as observed in the time resolved mappings, 
e.g., for T = 4K and hqH = 100 mT (Fig. [6) panels 
(iii) and (iv)). The step- like increase is followed by 
a slow relaxation characteristic of the time-scale t{ . 



9 




0.08 ± 0.01 s at T = 4 K to ra = 0.02 ± 0.005 s at 
T = 7.25 K. For iJ,oH = 100 niT, T3 decreases from 
Tg = 0.06 ± 0.01 s at T = 4 K to T3 = 0.02 ± 0.005 s 
at T = 6.75 K and for /io-ff = 135 mT, T3 de- 
creases from T3 = 0.045 ± 0.01s at T = 4K to 
T3 = 0.02 ± 0.005 s at T = 6.25 K. No signature of 
a crossover from the IMS to the Shubnikov phase 
is observed in the temperature and magnetic field 
dependence of t^. 

• The temperature dependence of T2 is character- 
ized by a linear decrease from T2 =0.047±0.01s 
at T = 4K to Ta = at T = 7.25K for a mag- 
netic field of fioH = 75 mT. For a magnetic field of 
fioH = 100 mT, T2 decreases from T2 =0.03±0.005 s 
at T = 4K to T2 = for temperatures above 
T = 6K. In contrast, for fioH = 135 mT, T2 = 
for all temperatures. Note, that r2 7^ in the IMS. 



V. INTERPRETATION 



FIG. 10: Panels (i) and (ii) show the values obtained for time 
constants ri , ri, T2 and rs as function of temperature T and 
magnetic field /j-qH. Note the different scaling for panel (i) 
and panel (ii). The lines serve as guide to the eye. 



Note, that the data points at low temperatures and 
low fields, where the kink is observed, are situated 
in the IMS. 

• In contrast, the characteristic shape of the relax- 
ation process exhibits a smooth increase with ex- 
ponential shape for high temperatures, namely for 
6.5K < T < 7.25 K at a magnetic field of fioH = 
75 mT, for 6 K < T < 6.75 K at ^o-ff = 100 mT and 
for 5 K < T < 6.25 K at ^qH = 135 mT. This is at- 
tributed to the lack of the above-mentioned sharp 
increase of intensity directly after the magnetic field 
is changed. 

We now concentrate on panels (ii) of Figs. [7||8] and |9] 
The decrease of intensity after the vertical line, labeled 
with trigger, is attributed to the VL being pulled out of 
the Bragg condition. The decreasing intensity is charac- 
terized by two different time-scales: T3 characterizes the 
overall time after the change of magnetic field direction 
when the scattering intensity has decreased to 1 /e of its 
initial value. T3 is indicated with a blue marker in each 
scan. In contrast, T2, indicated with a green marker, de- 
scribes the characteristic delay-time between the change 
of magnetic field direction and the response of the VL. 
The temperature and magnetic field dependence of T2 
and T3 is shown in Fig. 10 panel (ii): 



• T3 exhibits a linear decrease as a function of in- 
creasing temperature for all magnetic fields mea- 
sured: For fioH = 75 mT, T3 decreases from T3 = 



For an interpretation of our data, we first review the 
VL diffusion for uniform tilt before we qualitatively de- 
scribe the VL relaxation for a displacement of the mag- 
netic field direction. Following this we calculate the tem- 
perature and magnetic field dependence of the VL tilt 
modulus C44 using the VL diffusion model for uniform 
distortions. We finally consider the peculiarities of the 
relaxation process of the VL in the IMS. 



Vortex Lattice Elasticity and Diffusion for 
Uniform Distortions 



We have derived the characteristic relaxation rates Fi 
and r2 of VLs in section ([ll]) which are in the range of 
10~^ for Nb. In our experimental setup, the magnetic 
field direction was oscillated with a frequency of 0.2 Hz. 
Therefore, the associated relaxation process of the super- 
conducting VL can be calculated in the k = limit. For 
a change of the magnetic field direction, as used in our 
experimental setup, the relaxation process is essentially 
given by the VL tilt modulus C44(k = 0). 

For uniform distortions, the VL tilt modulus C44 = BH 
depends only on the applied magnetic field H which is in 
local equilibrium with the equilibrium induction . For 
the case that the spatially varying part of B is smaller 
than the average value of B, the highly nonlinear equa- 
tions of motion may be linearized. The response of the 
VL in a bulk sample (sample diameter ?■ >> A') to a 
changed magnetic direction field may then be written as 
damped diffusion process^^, using the diffusion equation 
derived by KeJ^. 

It is important to note that the VL initially responds 
to a change of applied magnetic field H solely at the 
surface of the superconductor, as the magnetic field is 
screened by supercurrents from the bulk of the specimen 
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(in particular by the Meissner domains in the IMS). Due 
to continuity conditions, the slope of the VL is slightly 
refracted at the surface of the sample. However, for the 
following description, the refraction of the magnetic field 
is neglected. 

According to the diffusion model, the distortion of the 
VL propagates from the surface of the sample into the 
bulk due to the finite elastic constants of the VL. The 
resulting diffusion equation of the tilt distortion u{x,t) 
of the VL is given by 



the resulting diffusion constant for a rotation of the mag- 
netic field with respect to the cylinder axis was calculated 
by BrandtP^. The resulting relaxation time is 



(2.405)2L) 



(19) 



where xq 
Jo{x). 



2.405 is the first node in the Bessel function 



du 
dt 



D 



d^ 
dx^ 



(13) 



with the diffusion constant D given by the ratio of the 
tilt modulus C44 and the viscosity 



D = 



C44 



(14) 



ry describes the viscous damping of the vortex motion by 
flux flow resistivity, assuming either vanishing pinning 
effects or an efficient depinning due to thermally assisted 
flux flow effects (cf. section (|ll]), eq. (10)). This yields 
for the diffusion constant D: 



D{T) 



(15) 



The diffusion equation can be solved easily, if the sam- 
ple is approximated by a flat plate with thickness 2r, 
ignoring the pre-existing field Hstat- We consider a con- 
ducting plate to which the transverse field i?osc is applied 
at t 0. The distribution of the field across the plate 
is then described through a square wave. As time pro- 
gresses, the edges and then the middle of the sample relax 
to the outside applied field. This corresponds to the so- 
called Dirichlet condition of the diffusion equation which 
is generally solved by a sum of cosine waves with a half 
period of 2r: 



'■ix,t) = 2^-Dn(cOS— ) 



n=l 

where n — 1,3,5... and 



2r 



2r 



riTTX 

fix) cos ax 



(16) 



(17) 



with the initial condition f{x,t = 0). 

The Fourier components thereby decay independently 
with a characteristic time oc 1/n^, so that after a short 
time, only the fundamental remains. This yields for the 
relaxation time for the mode with n ~ 1 



4r^ 



(18) 



For a sample of cylindrical shape and radius r with the 
magnetic field applied perpendicular to the cylinder axis. 



B. Vortex Lattice Relaxation Process 

We qualitatively describe the diffusion process of the 
VL. As the repetition cycle of -ffosc = 0.2 Hz is slow com- 
pared to the relaxation processes ti , t{ , T2 and ra of the 
VL, a complete relaxation can be presumed for each mea- 
surement cycle of the stroboscopic measurement, i.e. the 
VL is in the same equilibrium state before each change 
of magnetic field direction. 

We assume that the VL has relaxed for ^qHqsc — 
5 mT. The VL thus does not satisfy the Bragg condition. 
No scattering intensity is observed. We then consider 
the next cycle of the stroboscopic measurement where 
fJ-oHosc — 5 mT is decreased to /io^osc ~ mT: The VL 
firstly interacts with the changed magnetic field direction 
at the surface of the sample. The perturbation then dif- 
fuses into the sample. The VL relaxates into the Bragg 
condition. The relaxation of the VL from the equilib- 
rium position for /lo-ffosc = 5mT (and thus e = 2°) to 
fJ-oHosc = OmT (and thus e = 0°) yields a large displace- 
ment It 3> A' of the vortices. This relaxation process is 
characteristic of a slow exponential relaxation with the 
timescales ti and t{ as defined in the previous section. 
The temperature and magnetic field dependence of ti 



and r{ was given in Fig 10 panel (i 



We now assume that the VL has relaxed for hqHobc — 
mT, giving rise to maximum scattering intensity. If 
the direction of the magnetic field is shifted as iJosc is 
increased from /io^^osc =0 mT to 5 mT, the VL again in- 
teracts with the changed magnetic field direction at the 
surface of the sample. The perturbation diffuses into the 
sample, the VL is pulled out of the Bragg condition, and 
the scattering intensity at the detector decreases. Af- 
ter Ts, the intensity has decreased to 1/e of its initial 
value. T3 is thus a measure for the time-scale when the 
perturbation of the VL propagates across the complete 
sampl^^. 



C. Vortex Lattice Tilt Modulus C44 

According to the diffusion model fo r un iform distor- 
tions introduced by KePI and BrandPSl^ the VL tilt 
modulus C44 can be derived from ti and T3 using the 
eqs. (14), ( |Iol ), (15) and (19). The temperature depen- 
dence of Pn[T) is given by eq. (12). The equilibrium in- 



duction B is inferred from the modulus of the reciprocal 
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FIG. 11: Panel (i) depicts the temperature dependence of the 
equihbrium magnetization B for noH — 75 mT, 100 mT and 
135 mT. The hnes serve as guide to the eye. 



lattice vector [GvlI of the VL. Its temperature and mag- 
netic field dependence is given in Fig. 11 panel (ii). 

The temperature and field dependence of c1\{T,H), 
calculated from the measured values of ti, as well as 



cl^^{T,H) calculated from t[ and cl\{T,H), calculated 
from T3 is given in Fig. |12[ respectively. 

Only two data points could be obtained for C44(r, H) 
at noH = 100 mT for T = 4K and T = 6.5K from 
the time resolved mappings. C44 shows increasing vor- 
tex stiffness with increasing temperature, increasing from 
cl\ w 1- 10" TAm-i for T = 4K to cl\ w 1.4-104 TAm'^ 
forr = 6.5K. 

The VL tilt modulus c^^{T,H) shows increasing VL 
stiffness with increasing magnetic field. Moreover, for 

magnetic fields (i^H =75 mT and 135 mT, C44 exhibits a 
weak decrease by a factor of two for increasing tempera- 
ture from T = 4K to Tc. In contrast, for fioH =100 mT 

C44 = 0.3-10'* TAm"^^ shows no temperature dependence. 
For a temperature of T = 4K and a magnetic field of 

HqH =135 mT yields c^^l- 10'' TAm^i . For a tem- 
perature of T = 4K and magnetic fields of ijlqH =75 mT 

and 100 mT the VL tilt modulus c}^ « 0.3 - 10'' TAm"*. 

C44 (T, H) also shows increasing VL stiffness with in- 
creasing magnetic field. However, no significant temper- 
ature dependence is observed. Due to t( 3> T3, C44 yields 
cl\ « 5.5 - 10" TAm-i for ^i^H =135 mT, cj| « 2.5 - 10"* 
TAm-i for ^iqH =100 mT, and C4I « 2 - 10* TAj 



for 

HqH =75 mT. 

The expected temperature and magnetic field depen- 
dence of C44, calculated from literature values for Nb 
is given in Fig. 13 panel (i). The calculated diffusion 
time td is given in panel (ii). The effects of thermal 
depinning have been neglected. For a magnetic field of 
^qH =135 mT C44 assumes a value of 1.7-10* TAm~*, 
for hqH =100 mT C44 =0.8-10* TAm"* and finally for 
^iqH =75 mT C44 = 0.6 - lO^TAm"*. The temperature 
dependence of C44 is weak for all fields and reflects the 
temperature dependence of the equilibrium induction B. 

Qualitatively the measured values of cl\{T,H), 



■ / 01 1 1 1 








□ lOOmT 


1 ^44 


■ 135niT.- 


.---'J' 








C44 






-i — • — 



m 75mT 
□ lOOmT 
■ 135mT 



(ii) 



4 5 6 7 
T(K) 

Pulling VL out of Bragg condition 



" 6 
E 

^ 3 
o 























♦ M , ■ j 





■ 75mT 
□ lOOmT 

■ 135mT 



5 6 
T(K) 



FIG. 12: Panels (i) and (ii) show the calculated tempera- 
ture and magnetic field dependence of the VL tilt modulus 
clliT, H), cli{T, H) and c^KT, H), according to the model of 
BrandlP^. The lines serve as guide to the eye. 
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FIG. 13: The VL tilt modulus C44 calculated from literature 
values is plotted in panel (i), whereas panel (ii) yields the 
calculated timescale for the diffusion of a VL distortion in a 
cylindric Nb sample, according to the model of BrandtP^. The 
lines serve as guide to the eye. 
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cl\{T,H) and c^\{T,H) compare well with the theoret- 
ical value of C44{T,H): Firstly, the magnetic field be- 
haviour is consistent with an increased VL stiffness C44 
for increasing field, according to eq. ([9|. The most ac- 
curate agreement is obtained for cl\{T,H). The origin 
of the deviation is most likely due to uncertainties of the 
extrapolation of the normal conducting resistivity pn to 
low temperatures, associated with the extrapolation of 

the RRR. The values of cl\{T,H) are lower by a factor 
two from the expected values. In contrast, the values for 
cl\{T, H) exceed the calculated values by a factor of 3.5. 
This deviation is presumably resulted by the definition 
used for t[ and t^. 

We have introduced in the previous paragraphs that 
the Fourier modes of the VL relaxation decay indepen- 
dently with a time constants cx 1/n^. Note, that the 
values of ti exceed T3 by a factor of approximately two 
to three. Whereas the fundamental relaxation mode ri 
may be identified with the Fourier mode for n = l^^Si, the 
values obtained for T3 suggest that the corresponding re- 
laxation process is a mixture of several different mech- 
anisms. Moreover, the involvement of the shear mod- 
ulus C66 in the relaxation process is still unclear. The 
shear modulus is responsible for the perfection of the lo- 
cal structure of the VL. However, no azimuthal smearing 
of the scattering pattern was observed during the relax- 
ation process of the VL, indicating a disordering akin to 
a melting transition of the VL. This leads to the follow- 
ing picture of the VL relaxation, already accounting for 
most the observed features: 

(i) In general, a faster propagation of perturbations 
in the VL is observed for increasing temperatures and 
increasing magnetic fields. This general behaviour is ex- 
plained by the decreased damping r] for increasing tem- 
perature, according to an increase of Pn{T) with T as 
given in eq. (fTol) and increasing VL stiffness C44 with 



increasing field according to C44 = BH. 

(ii) The time resolved mappings show that at high 
temperatures and high magnetic fields in the Shubnikov 
phase, the VL responds to the changed magnetic field 
direction as a single, stiff lattice with its motion charac- 
terized by the diffusion constant D. A single, exponential 
relaxation is thus observed. 

(iii) In contrast, due to the low stiffness and strong 
damping for low temperatures and fields in the IMS, the 
macroscopic relaxation is typical of a slow exponential re- 
laxation with large mosaic spread which is characterized 
by an additional fast process on a short timescale. This 
is attributed to the decomposition of the VL into Shub- 
nikov domains and Meissner phase domains in the IMS. 
The origin of the additional fast process will be discussed 
below. 

(iv) The macroscopic relaxation of the VL between the 
two equilibrium positions is associated with a large vor- 
tex displacement. It is strongly dependent on the VL 
topology, as the crossover from the IMS to the Shubnikov 
phase is reflected in the mere temperature dependence of 



(v) Ta is a measure for the time of a VL distortion 
propagating through the sample. No signature of the 
transition from IMS to the Shubnikov phase is observed 
in the temperature and magnetic field dependence of T3 
which is thus insensitive to the VL topology. 



D. Vortex Lattice Relaxation in the Intermediate 
Mixed State 

Certain characteristic features show up exclusively for 
data points in the IMS. Typical data of the relaxation 
process for T — AK and ij^qH =100 mT is shown in 



Fig. 14 The time range displayed corresponds to a short 
time-scale after the magnetic field direction was changed. 
In the time resolved mapping (panel (i)), a sharp increase 
of mosaicity shows up which is associated with a sharp 
increase of intensity observed in the fixed angle scans for 
the equilibrium position at 4> = 2.75° (panel (ii)). We 
note, that the characteristic time-scale of this feature is 
well below for all temperatures and fields in the IMS. 

We have introduced T3 as a measure for the time, nec- 
essary for a perturbation of the VL to cross the sample. 
This locates the related relaxation process responsible 
for the sharp increase of intensity and mosaicity at the 
surface of the sample. This effect is hence attributed to 
branching of the Shubnikov domains at the surface of the 
sample in the IMS. The branching of Shubnikov domains 
leads to a fine VL structure at the surface, consisting 
of connected Shubnikov domains with open topology en- 
closing Meissner islands. In particular, no rigid VL is ob- 
served. Branching of the VL is responsible for the large 
intrinsic mosaicity as well. 

As remaining feature, we now discuss the character- 
istic delay T2 when the VL is pulled out of the Bragg 
condition. A delay T2 =30 ms is visible in both the time 
resolved mapping (panel (i)) and the corresponding fixed 
angle scan for = 0.5° (panel (ii) of Fig. 14). Note, that 
T2 7^ solely for data points in the IMS. A sharp increase 
of intensity is observed simultaneously at the new equilib- 
rium position (0 = 2.75°) for these measurement points. 
This sharp increase was attributed to a mechanism close 
to the surface of the sample. The IMS is characterized 
a bending of vortices at the surface. The scattered in- 
tensity exactly at the position q = Gvl = ^ is 
caused by the VL which points to the initial equilibrium 
position, thus which is not bent due to branching. T2 
thus describes the delayed reaction of the unbranched 
VL buried beneath the surface of the sample. 



VI. CONCLUSIONS AND OUTLOOK 

In our study we measure the dynamic properties of 
bulk VLs with time resolved stroboscopic SANS com- 
bined with a time varying magnetic field setup. We show 
that the qualitative magnetic field and temperature de- 
pendent behaviour of the VL in the Shubnikov phase for 
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FIG. 14: Panel (i): time resolved mapping for T =4K and 
fioH =100 mT. The relaxation process of the VL is shown 
for a short time-scale after the magnetic field direction was 
switched. Panel (ii) shows the associated fixed angle scans 
for <j) — 0.5° and (p — 0.2.75° as indicated with the vertical 
broken red lines in panel (i). 



uniform distortions can be described in reasonable agree- 
ment with a diffusion model by Brandt^ and KeJ^. 
The values obtained for the VL tilt modulus C44, diffu- 
sivity D and damping r/ can be reproduced theoretically. 
We further argue that the topology of the VL is reflected 
sensitively in the associated diffusion process. This is 
readily seen for measurement points in the IMS where a 
second, very fast relaxation process maybe attributed to 
branching of the VL domains at the surface of the sam- 
ple. Such a branching was observed by previous high res- 
olution decoration techniques on bulk Nb single crystal 
sample^^^l^Sll xhe effect of surface treatment and rough- 



ness on the details of the Landau branching remains as 
a subject of future studies. 

Our study represents a experimental technique for how 
to access directly VL melting and the formation of vortex 
glass states in unconventional superconductors, notably 
the cuprates, heavy-fermion systems, borocarbide or iron 
arsenide systems. The possibility to precisely determine 
the pinning properties of vortices in future experiments 
on samples of varying purity is of great relevance for the 
research on technical applications of superconducting de- 
vices, as the pinning properties are intimately related to 
the maximal critical current density of superconductors. 

Furthermore, the technique developed for our study 
is of general relevance for materials exhibiting complex 
forms of magnetic order, i.e. long range helical or der as 
observed in materials without inversion symmetrjHSm]^ 
magnetic Skyrmion lattices, as observed recently in the 
helimagnet MnSp^lor colloidal magnetic suspensions, liq- 
uid crystals and also Bose condensates and glasses of 
magnetic triplet excitations in quantum magnet^^. 

To further increase the time resolution, the TISANE 
technique can be used instead of stroboscopic small angle 
neutron scatterin^^. TISANE benefits from a neutron 
chopper which is placed upstream of the sample position 
at the distance Li to the sample. By carefully adjust- 
ing the distances L2 and Li, the chopper and control 
parameter duty cycle frequencies and phase, a coherent 
summation of the scattered neutrons at the detector po- 
sition in space and time without time smearing can be 
achieved. The accessible time resolution of TISANE is 
mainly determined by the time resolution of the detector 
and the opening time of the chopper system. Timescales 
from ^s to several ms are possible, closing the gap to the 
inelast ic tech nique NRSE. For a detailed description, we 
refer tcPHl. 
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